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The third law of thermodynamics constrains the phase diagram of systems with a first-order 
quantum phase transition. For zero conjugate field, the coexistence curve has an infinite slope at 
T = 0. If a tricritical point exists at T > 0, then the associated tricritical wings are perpendicular 
to the T = 0 plane, but not to the zero-field plane. These results are based on the third law 
and basic thermodynamics only, and are completely general. As an explicit example we consider 
the ferromagnetic quantum phase transition in clean metals, where a first-order quantum phase 
transition is commonly observed. 
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Systems with 


First-order phase transitions are ubiquitous in nature, 
the solid-to-liquid and liquid-to-gas transitions being the 
most commonly observed ones. Another common exam¬ 
ple of a first-order transition is the ferromagnetic tran¬ 
sition below the Curie temperature as a function of an 
external magnetic field. First-order transitions are char¬ 
acterized by a coexistence curve in the phase diagram 
along which the two phases coexist in thermodynamic 
equilibrium. (The coexistence curve may be the projec¬ 
tion of a higher-dimensional coexistence manifold into a 
particular plane in the phase diagram.) 

It has long been known that the curvature of the coex¬ 
istence curve is determined by the discontinuities of cer¬ 
tain observables across it. The Clapeyron-Clausius (CC) 
equation relates the slope of the coexistence curve in the 
pressure-temperature ( p-T ) plane to the discontinuities 
of the entropy and the volume jT]: 


/ dp\ As 
\dT) H A v 


(1) 


where As = Si — S 2 and Av = V\ — V 2 with Sq 2 and 
the specific entropy and volume per particle, respectively, 
in the two phases. For definiteness, let 1 and 2 label the 
ordered and disordered phases, respectively, and for later 
reference we indicate that an appropriate external field 
H, if any, is held constant in taking the derivative. 

The CC equation (|T|) and its analogs in different planes 
of the phase diagram are very general, as they rely only 
on basic thermodynamic arguments. In this Letter we 
show that for quantum phase transitions, when combined 
with the third law of thermodynamics, they provide in¬ 
teresting constraints on the shape of the phase diagram. 
We will consider a pressure-driven transition at T = 0 
that is first order, remains first order at low T, and turns 
second order at higher T via a tricritical point (TCP). 



Figure 1: Schematic phase diagram showing a line of first- 
order transitions at low T separated from a line of second- 
order transitions at higher T by a tricritical point (TCP). 
In a nonzero conjugate field H tricritical wings emerge from 
the TCP. These are surfaces of first-order transitions that are 
bounded by lines of second-order transitions and terminate 
in two quantum wing-critical points (QWCP) in the T = 0 
plane. 


The schematic phase diagram in the space spanned by 
T, p, and H, where H is the field conjugate to the or¬ 
der parameter, is shown in Fig. |T] As we will see, the 
detailed shape of this phase diagram at low T is con¬ 
strained by thermodynamics. Our arguments leading to 
this conclusion are completely general; however, as an 
explicit example we will discuss the quantum ferromag¬ 
netic transition in clean metals |2|. Another example of a 
first-order quantum phase transition with a TCP in the 
phase diagram is the Ising antiferromagnet dysprosium 
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aluminum garnet fJ]. 

We are interested in a system with T, p, and H as in¬ 
dependent variables. Denoting the order parameter by 
M , the appropriate thermodynamic potential is the gen¬ 
eralized Gibbs free energy |3J 

G = U -TS+pV - HM 

= 7 iN , (2a) 

whose differential is 

dG = —SdT + Vdp - MdH + gdN . (2b) 

Here S, V, and fi are the system’s entropy, volume, and 
chemical potential, respectively, and N is the particle 
number. From Eqs. ([2]) we obtain the Gibbs-Duhem re¬ 
lation 


dg = dg = —sdT + vdp — mdH 


( 3 ) 


where g , s, v, and m are the generalized Gibbs free en¬ 
ergy, entropy, volume, and order parameter per particle, 
respectively. On the coexistence curve the chemical po¬ 
tentials of the two phases must coincide. Using this con¬ 
dition with Eq. |3} at fixed external field leads to Eq. ([I}. 
An analogous argument yields 


dT\ 

dH)> 


Am 

As 


( 4 ) 


The CC equations are completely general. When 
applied to a quantum phase transition, the third law pro¬ 
vides the additional constraint A s(T —> 0) —► 0. To be 
specific, let us assume that in either phase the entropy 
vanishes as s(T —> 0) = 7 T". In particular, if the phases 
are Fermi liquids (see below) then n = 1 and 7 is the 
specific-heat coefficient. For asymptotically low temper¬ 
atures we thus have 


c£T N 

d P / h 
dT_\ 
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(5a) 

(5b) 


In addition, we obtain from the equilibrium condition in 


conjunction with Eq. (2b I a third CC equation 


dH \ Av 


dp J r 


Am 


(5c) 


These three CC equations are the basis of our discussion. 

Let us start by briefly discussing the obvious coexis¬ 
tence region in the T-p plane, which is labeled “ordered” 
in Fig. [T] The first-order transition across this plane, 
which is driven by the external field, does not involve 
any change in either entropy or specific volume. We 


thus have As = Av = 0. Equations (5b 5c) then imply 


( dH/dT) p = ( dH/dp)T = 0. This identifies the H = 0 


plane as the locus of the coexistence curves. (dT/dp) h= 0 , 
Eq. (5aI, is indeterminate, which is consistent with the 
fact that any curve in the H = 0 plane below the tran¬ 
sition temperature is a coexistence curve. The CC equa¬ 
tions thus correctly describe the coexistence plane, but 
do not provide any nontrivial information. 

This changes as we consider the other coexistence sur¬ 
faces, viz., the tricritical wings. Obviously, we have 
Am > 0 across any first-order transition, but now As 
and Av will not be zero. To find Av, we turn to scaling 
theory. Scaling is often thought of as valid only at second- 
order transitions. However, Fisher and Berber [5| have 
shown that finite-size scaling considerations allow for the 
definition of a diverging length scale even at a first-order 
transition. Consequently, a classical first-order transition 
can be considered a limiting case of a second-order transi¬ 
tion, and the homogeneity laws, exponent relations, etc., 
that are known from the scaling description of second- 
order transitions still hold. This formalism has recently 
been generalized to the case of quantum first-order tran¬ 
sitions [ 6 ], and we now apply it to the problem under 
consideration. Let r = (p — p*)/p* be the dimensionless 
distance from the transition at T = 0. Then the general¬ 
ized Gibbs free energy obeys a homogeneity law [ 6 j 

g{r, H, T) = b~ (d+z) (r b d+z , H b d+z ,Tb z ) . ( 6 ) 


Here b > 0 is the length scaling parameter, <f>g is a 
scaling function, and we have made use of several expo¬ 
nent values that characterize a first-order quantum phase 
transition (see Ref. [6] for details): z is the relevant dy¬ 
namical exponent | 7 j, the inverse correlation length ex¬ 
ponent has its largest possible value \/v = d+z, and the 
scale dimension of the field, [H] = fi5/v, reflects the fact 
that the order-parameter exponents are /3 = 0, 8 = 00 , 
such that /3S = 1. This is a generalization of the scal¬ 
ing description of classic first-order transitions given by 
Fisher and Berker (5,. Differentiating g with respect to 
r oc p, we see that the scaling part of the specific volume, 
v = dg/dp = dg/dr, has a zero scale dimension. This 
implies a discontinuity of the specific volume across the 
transition, and a corresponding ^-function contribution 
to the compressibility k = —(dv/dp)/v. This is in direct 
analogy to the latent heat at a classical first-order tran¬ 
sition and can be interpreted as a “latent volume” at a 
pressure-driven QPT, i.e., the volume changes by a finite 
amount upon an infinitesimal change in pressure [B]. We 
also note that differentiating with respect to r again we 
obtain k oc |r| _1 ; see Ref. |5] for an interpretation of this 
power-law divergence at a first-order transition in terms 
of finite-size scaling. The compressibility is positive in 
thermodynamic equilibrium, and we thus have Av > 0 
across any coexistence curve contained in the wings. The 
only remaining question is the sign of Ay. For a transi¬ 
tion from an ordered phase to a disordered one, one would 
naively expect Ay < 0. We will adopt this expectation 
for now and give a more detailed discussion below. 
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Figure 2: Measured phase diagram of UGe 2 , with the same 
notation as in Fig. [T] The main figure is adapted from Ref. llOt 
the squares represent data points, the lines and surfaces are 
guides to the eye. The inset shows the H = 0 plane with data 
from Ref. flU Note the extremely steep drop of the Curie tem¬ 
perature past the tricritical point. The ferromagnetic phase 
consists of two phases, FM1 and FM2, separated by a line 
of first-order transitions at low temperatures that ends in a 
critical point (CP). 


From Eq. (5a I we conclude that ( dT/dp)n < 0, i.e., 
the slope of any coexistence curve at fixed H is negative, 
and it approaches — oo as 1/T" for T —» 0. Similarly, 
Eq. (5b) yields ( 8T/dH) p > 0, and it approaches +oo as 
1/T n for T —> 0. This means that the wings are necessar¬ 
ily perpendicular to the T = 0 plane, and, in particular, 
the coexistence curve in zero field has an infinite slope at 
the QPT. Finally, Eq. (5c I implies that ( dH/dp)T is pos¬ 
itive and finite, which implies that the wings are tilted 
in the direction of the disordered phase and are not per¬ 
pendicular to the p-axis. 


We now turn to an explicit example that illustrates all 
of the above considerations, namely, the quantum phase 
transition in clean metallic ferromagnets mm- In this 
case, the exponent n in Eqs. [5]is n = 1, and the dynam¬ 
ical exponent z in Eq. (JgJ) is 2 = 1. There is a second 
dynamical exponent z = 3, but for our purposes z = 1 
yields the dominant contribution (see Refs. El m for a 
detailed discussion of this point). The order parameter 
to is the magnetization, H is the external magnetic field, 
and the phase diagram is generically observed to have the 
topology shown in Fig. [I] The features discussed above 
are indeed universally observed in all cases where the tri¬ 
critical wings have been mapped out in detail. As an 
example, we show the experimentally determined wings 
in UGe 2 in Fig. [2j for other examples, see Ref. EJ Note 
the extremely sharp drop of the Curie temperature for 
pressures above the tricritical pressure that is apparent 
in the inset. Also of interest is the transition inside the 
ferromagnetic phase (from FM1 to FM2), which is first 


order at low temperatures. Our considerations apply to 
this transition as well, and the steep drop of the transi¬ 
tion temperature is again consistent with an infinite slope 
of the coexistence curve at T = 0. 

We now return to the issue of the sign of Ay, or more 
generally As, across the coexistence curve. From Eq. ^ 
we see that ( ds/dp)T,H = ~(dv/dT) Pt H = ~va p , with 
a p = ( dv/dT) p /v the thermal expansion coefficient. An 
increase in entropy with increasing pressure thus implies 


a p < 0. Returning to Eq. (5aI, and remembering that 


Av > 0 since the compressibility is necessarily positive, 
we see that a decreasing Curie temperature with increas¬ 
ing pressure implies a negative thermal expansion coeffi¬ 
cient, and vice versa. Consistent with this, the thermal 
expansion coefficient at low T is indeed negative in UGe 2 
[TE] , MnSi |T3j > and ZrZn 2 [3] , which all are low-T ferro¬ 
magnets with qualitatively identical phase diagrams. It 
is interesting that a p < 0 by itself implies that the high- 
pressure phase must be the paramagnetic one. We also 
note that the volume is discontinuous if the coexistence 
curve is crossed at fixed p as a function of T as well as 
at fixed T as a function of p. This is intuitively obvi¬ 
ous and also follows from Eqs. Accordingly, a p 

at the first-order transition has a (5-function contribution 
that reflects the same “latent volume” as the correspond¬ 
ing (5-function contribution to the compressibility. This is 
consistent with the experiment by Kabeya et al PE], who 
observed a pronounced negative peak in a p at the tran¬ 
sition, which they attributed to a broadened first-order 
transition. 

While in most quantum ferromagnets hydrostatic pres¬ 
sure destroys the ferromagnetic order, there are systems 
in which the opposite occurs. An example is YbCu 2 Si 2 , 
which is paramagnetic at ambient pressure, but becomes 
ferromagnetic upon the application of hydrostatic pres¬ 
sure of roughly 10 GPa m , see Fig. [3] The thermody¬ 
namic arguments presented above then predict that a p 
in this material, at this pressure and at low temperature, 
must be positive. The thermal expansion data of Ref. [TB] 
are consistent with this prediction, although not quite 
conclusive, as they focused on a higher temperature re¬ 
gion. There also are materials where hydrostatic pressure 
drives the system away from ferromagnetic order, while 
uniaxial stress favors it; an example is UCoAl [T7]. This 
can be understood by realizing that in many solids the 
thermal expansion coefficient is anisotropic to the point 
of being positive along some crystal axes, but negative 
along others }T5], 

We finally briefly discuss an explicit equation of state 
that has been used to describe the qualitative phase di¬ 
agram of metallic quantum ferromagnets and leads to a 
schematic phase diagram as shown in Fig. 0C2] . It is 
derived by minimizing a generalized Landau functional 


/(to) = —Hm + rrn 2 + wm 4 ln(m 2 + T 2 ) + ' 


(7) 


with respect to the magnetization to. Here to, H , and 
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Figure 3: Temperature-pressure phase diagram of YbCu 2 Si 2 - 
Ferromagnetism is induced by hydrostatic pressure p > P c ~ 
8.25 Gpa. From Ref. 1151 



T are measured in suitable microscopic units, and r, u, 
and w are parameters of the generalized Landau theory. 
The physical origin of the term with coupling constant 
w is due to soft fermionic excitations that couple to the 
magnetization; this has been discussed in detail before [S] 
J5DJ12I] and will not be repeated here. In zero held, H = 0, 
the logarithmic term leads to a first-order transition at 
r = r\ = w e~ x ~ u ! w where the magnetization value is 
mi = yfr\, and to a tricritical point at T tc = e~ u ^ 2w . The 
coexistence curve can easily be obtained explicitly DU; 
here we just quote the asymptotic behavior for r —> r±, 

T(r —*■ tt) = -4= (tt — r) 1 / 2 (8) 

Vw 

which yields 

dT/dr\ r ^ ri = -1/2 wT , (9) 


in agreement with the general conclusions drawn above 


from thermodynamics, see Eq. (5a I and the related dis¬ 


cussion. In the presence of a small magnetic held, one 
finds for the coexistence curve in the T = 0 plane 

H = mi (1 + \ 5r + 0(Sr 2 ) , (10) 

\ 11 w J 

where Sr = r — r\. This reflects the linear slope of the 
tricritical wings with respect to the r-axis that follows 


from Eq. (5cI . A more involved, but elementary analysis 
shows that the tricritical wings are perpendicular to the 
T = 0 plane everywhere. We stress that these proper¬ 
ties are not tied to the specihc physical mechanism that 
underlies the free-energy function ([Tj; they must be true 
for any model that leads to a hrst-order quantum phase 
transition and correctly rehects thermodynamics. 

We conclude with some additional discussion points. 
(1) The most often observed shape of the phase diagram 


in quantum ferromagnets, with increasing hydrostatic 
pressure driving the system into the disordered phase, 
is not what one might naively expect. In a fluid analogy, 
this is equivalent to what is observed in H 2 0 and H 2 S, 
while in most fluids increased pressure stabilizes the or¬ 
dered phase. In quantum ferromagnets the latter can 
also occur, see Fig. [3j but it is not common. As we have 
shown, this feature of the phase diagram is tied to the 
sign of the thermal expansion coefficient, which tends to 
be negative in low-temperature ferromagnets. Regardless 
of whether hydrostatic pressure induces or destroys fer¬ 
romagnetism, Eq. (5c) implies that the tricritical wings 
must always extend in the direction of the paramagnetic 
phase. (2) All actually measured tricritical wings show all 
of the structural features discussed above, as they must, 
since the former hinge on basic thermodynamics only. 
However, occasionally schematic drawings of wings in 
the literature violate these thermodynamic requirements: 
They show wings that are perpendicular to the pressure 
axis and/or not perpendicular to the T = 0 plane. (3) 
The discontinuity of the specific-heat coefficient across 
the first-order transition, 0 < Ay < oo, has an inter¬ 
esting implication for the dynamical critical exponents 
in the system. Since the static and dynamic specific-heat 
exponents a and z c are related by a hyperscaling relation 
a = v(z c — d), with v the correlation-length exponent and 
d the spatial dimensionality of the system |B], a discon¬ 
tinuous specific-heat coefficient (a = 0) implies z c = d. 
This is indeed the value of z c within an explicit theory 
that describes the first-order transition, see Ref.® 
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